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ABSTRACT 

A neutron star in a long-Uved, low-mass binary can easily accrete enough matter to 
replace its entire crust. Previous authors noted that an accreted crust, being formed 

^-^ . from the burning of accreted hydrogen and helium, allows a series of non-equilibrium 

reactions, at densities ^ 6 x lO^^gcm"^, which release a substantial amount of heat 

Qv I (~ IMeV per accreted nucleon). Recent calculations by Schatz et al. showed that the 



crystalline lattice of an accreted crust is also likely to be quite impure. This paper 



CIh' discusses the thermal structure of such a neutron star and surveys how the crust re- 

Q ■ actions and impurities affect the crust temperature. During accretion rapid enough to 

^— > ■ make the accreted hydrogen and helium burn stably (M ~ lO^^AfQyr^^; typical of 

C^ . the brightest low- mass neutron star binaries), most of the heat released in the crust is 

conducted into the core, where neutrino emission regulates the temperature. As a result 
k>>( , there is an inversion of the thermal gradient: the temperature decreases with depth in 

H ' the inner crust. The thermal structure in the crust at these high accretion rates is 

insensitive to the temperature in the hydrogen/helium burning shell. When the crust 
is very impure, the temperature can reach ~ 8 x 10^ K at densities ^ 6 x 10^^ gcm~^. 
This peak temperature depends mostly on the amount of heat released and the thermal 
conductivity and in particular is roughly independent of the core temperature. The 
high crust temperatures are sufficient to melt the crystalline lattice in thin layers where 
electron captures have substantially reduced the nuclear charge. 

Subject headings: accretion, accretion disks — conduction — stars: neutron — stars: 
interiors 
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1. Introduction 

Many of the known neutron stars reside in low-mass x-ray binaries. These sources typically 
accrete at rates of 10^^"*^ Mq yr^"*^ to IO^^Mq yr~^ and show no conclusive evidence, such as cyclotron 
lines or coherent pulsations (in the persistent emission), of a magnetic field. In contrast to studies of 
isolated neutron star cooling (for a review, see Tsuruta 1998), there has been much less interest in 
the interior thermal state of an accreting neutron star. Originally, studies of how accretion affected 
a neutron star's thermal structures were motivated by the challenge of explaining the type I x-ray 
bursts of some of these sources. Both Lamb &: Lamb (1978) and Ayasli & Joss (1982) estimated the 
steady-state core temperature by balancing the heating from hydrogen and helium shell burning 
with neutrino and radiative losses. Later, Fujimoto et al. (1984) and Hanawa & Fujimoto (1984) 
calculated the thermal evolution of the entire neutron star, both for steady hydrogen and helium 
burning (Fujimoto et al. 1984) and for repeated shell flashes (Hanawa & Fujimoto 1984). In all of 
these works, the only heat sources considered were the hydrogen/helium burning and the influx of 
entropy by the accreted matter. Both Fujimoto et al. (1984) and Hanawa & Fujimoto (1984), who 
considered accretion rates ^ 5 x 10~^MQyr~^, found that the deep crust and core would gradually 
become isothermal, at a temperature > 10® K, if there were no enhanced neutrino cooling; otherwise, 
the core would remain chilled at temperatures ~ 10 K. Without heat sources in the crust, the 
temperature of the deep crust tracks that of the core and is therefore sensitive to the cooling from 
neutrino processes active in the core. 

Unlike an isolated neutron star, the crust of an accreting neutron star is not in statistical 
nuclear equilibrium, but rather has a composition set by the nuclear history of the accreted material 
(Sato 1979; Blaes et al. 1990, 1992; Haensel & Zdunik 1990a, b). The atmosphere is composed of 
the accreted helium and hydrogen and any metals present (see Bildsten, Salpeter, & Wasserman 
1992 for a discussion). This accumulated fuel eventually burns to heavier elements. The accretion 
of fresh fuel shoves the original crust deeper and, if continued over a long enough interval, will 
eventually replace the original crust with one formed by the ashes of hydrogen/helium burning. 
Compression of the crust by the weight of continually accreting material induces non-equilibrium 
reactions that release heat. 

The composition of the replaced crust is uncertain. Improved treatments of the physics of 
hydrogen and helium nuclear burning revealed that the ashes of this burning are unlikely to be a pure 
species, e.g., iron. The mixture is formed by the rp-process (Wallace &: Woosley 1981; Champagne 
& Wiescher 1992; Van Wormer et al. 1994; Schatz et al. 1998), a sequence of rapid proton captures 
onto seed nuclei provided by helium and CNO burning. Calculations of the nucleosynthetic yield 
from the rp-process have been done, both for unstable burning during an x-ray burst (Koike et al. 
1999) and for steady-state hydrogen and helium burning (Schatz et al. 1999). The ashes of the 
stable burning are a motley mix of iron-peak elements, so that the crust formed from these ashes 
will likely be very impure (Schatz et al. 1999). 

This paper studies the crustal temperatures of steadily accreting neutron stars with low mag- 



netic fields. There are two differences witli Miralda-Escude, Paczynski, & Haensel (1990) and 
Zdunik et al. (1992), botli of wliich included crust reactions in the neutron star's thermal balance. 
First, this work considers the stable regime of hydrogen/helium burning, which requires rapid accre- 
tion (near the Eddington limit, ~ lO~*M0yr~^). At the low accretion rates (M < lO~^''M0yr~^, 
roughly two orders of magnitude less than the Eddington limit) considered by Miralda-Escude et al. 
(1990) and Zdunik et al. (1992), the crust is basically isothermal, with a temperature locked to 
that of the core. Second, this work allows for an impure crust by surveying both high- and low- 
conductivity cases. Previous calculations have assumed the impurity concentration to be much less 
than unity. 

At the rapid accretion rates considered here, the neutrino luminosity from modified Urea 
processes and crust bremsstrahlung is significant and causes the temperature to decrease with 
depth in the inner crust (Brown & Bildsten 1998). Almost all of the heat produced in the crust 
flows inward. Moreover, the reduced conductivity of the impure crust produces a peaked thermal 
profile with a maximum temperature where the nuclear reactions heat the crust, at densities > 
6 X lO^^gcm^^. The thermal profile in the crust is primarily determined by the ability of the 
inner crust to conduct a flux of ~ 1 MeV per accreted baryon into the core, and is insensitive 
to the temperature of the hydrogen/helium burning shells. If the crust is very impure, the crust 
reaches temperatures ~ 8 x 10* K; the value of this temperature only weakly depends on the core 
temperature. 

Electron captures in the crust reduce the charge of the nuclei (Z) and hence the electrostatic 
binding of the lattice. For the hottest temperatures in the crust, this low-Z lattice melts where the 
charge is lowest (Z < 15 for the composition of Haensel & Zdunik 1990a). As a result, the inner 
crust of the neutron star comes to resemble a "layer cake," with alternating layers of lattice and 
liquid. 

This paper is relevant for the brightest low-mass x-ray binaries. These weakly magnetized 
neutron stars are considered possible progenitors of millisecond pulsars (for a review, see Bhat- 
tacharya 1995), and there has been much theoretical interest in the evolution of the crust magnetic 
field (Romani 1990; Geppert & Urpin 1994; Urpin & Geppert 1995; Konar & Bhattacharya 1997; 
Brown & Bildsten 1998; Urpin, Geppert, &: Konenkov 1998). Many of these neutron stars rotate 
within an apparently narrow range of spin frequencies ~ 300 Hz (van der Klis 1998, e.g.). One 
possibility for this convergence of spin frequencies is that gravitational radiation from the neutron 
star balances the accretion torque (Bildsten 1998; Andersson et al. 1999). The source for the grav- 
itational radiation could be a mass quadrupole formed by misaligned electron capture layers in the 
crust (Bildsten 1998) or a current quadrupole from an r-mode instability in the core (Andersson 
1998; Friedman & Morsink 1998; Lindblom, Owen, & Morsink 1998; Owen et al. 1998; Andersson, 
Kokkotas, &: Schutz 1999; Andersson, Kokkotas, & Stergioulas 1999). All of these problems depend 
on the thermodynamics of the neutron star's crust and core and motivate this paper. 



1.1. An overview of the problem 

A neutron star has several distinctive regions. The core consists of uniform npe^ (in the least 
dense parts). At a baryon density less than n ~ CGn^, where n^ = 0.16 fm^ is the saturation 
density^ of nuclear matter, individual nuclei appear (Pethick, Ravenhall, & Lorenz 1995). The 
portion of the neutron star exterior to this point, the inner crust, is composed of nuclei, degenerate 
neutrons, and relativistic degenerate electrons. Where the electron Fermi energy is less than about 
twice the nuclear bulk energy (~ 30MeV; see Pethick &: Ravenhall 1998), free neutrons can no 
longer exist in /^-equilibrium. This point, neutron drip, has a density ~ 0.0023ns (4 x 10~^ fm~ ) 
and marks the boundary between the inner and the outer crust. At lesser densities, the crust is 
made of nuclei and electrons, with the degenerate electrons supplying the pressure. The outer and 
inner crust collectively occupy the outermost kilometer or so of the neutron star and contain a total 
mass of order O.OIMq. 

The boundaries of the crust are demarcated by surfaces of constant pressure. For a thin 
crust, the mass above a given isobar is fixed by the surface gravity and area. As a consequence, 
accretion during a time brief compared with M/M (so that the overall structure of the neutron 
star remains roughly constant) pushes the underlying crust through these compositional boundaries. 
Low-mass x-ray binaries live for more than 10^ yr (Webbink, Rappaport, & Savonije 1983); even 
at M = lO^^'^MQyr"-'^, the neutron star can easily accrete enough material from the secondary 
to replace its entire crust. This replaced crust is composed of the ashes of hydrogen and helium 
burning and is quite different in composition from the original. Unlike during the neutron star's hot 
birth, the crust does not burn to nuclear statistical equilibrium. As the ashes of the hydrogen and 
helium burning are pushed deeper into the crust, the rising electron Fermi energy induces a series of 
electron captures (Haensel & Zdunik 1990b; Sato 1979; Blaes et al. 1990). Further compression of 
this low-Z material causes neutron emissions and pycnonuclear reactions (Haensel & Zdunik 1990b; 
Sato 1979). A schematic of the composition, from the calculation of Haensel & Zdunik (1990b), 
is shown in Figure 1. Each decrease in the nuclear charge Z (bottom panel) is from an electron 
capture, and each decrease in the nuclear mass number A (top panel) is from a neutron emission. 
Where a pycnonuclear reaction occurs, both Z and A double. 

The pycnonuclear reactions and neutron emissions liberate E'at ~ 1 MeV for each baryon 
accreted and heat the crust at a rate L]\f ~ lO^^(M/lO~^M0yr~"^) ergs""*^. To set the scale for the 
core temperature, note that the core must be at a temperature ~ 4 x 10^(-Ljv/10^^ ergs^^)^'*^ K 
for modified Urea processes (Friman &: Maxwell 1979; Yakovlev & Levenfish 1995) to radiate a 
neutrino luminosity equal to Ljy. If the nucleons in the core were superfluid and the modified Urea 
processes correspondingly suppressed, crust neutrino bremsstrahlung (Maxwell 1979; Yakovlev & 
Kaminker 1996; Haensel, Kaminker, & Yakovlev 1996; Itoh et al. 1996) can also balance Li\i for an 



^Density is measured in units of fm '* and pressure in units of MeVfm '^. Note that throughout the crust, 
p ~ nniu ~ 1.66 x 10^''(n/0.1 frn^'^) gcm^'\ where m„ = 1.66 x 10^^'* g is the atomic mass unit, and 1 MeVfm^'* — 
1.6 X lO^^dynecm"^ 
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Fig. 1. — Composition of an accreted crust, as computed by Haensel &: Zdunik (1990b). Shown 
are the mass number (top panel) and charge number (bottom panel) as a function of pressure. The 
mass resting above a given isobar is given by equation (8). Electron captures produce a downward 
jump in Z; neutron emissions produce a downward jump in A. The upward jumps in both Z and 
A indicate a pycnonuclear reaction. 




inner crust temperature ~ 6 x 10^(LAr/10^^ ergs ^)^/^ K. 

The luminosity from the steady-state hydrogen/hehum burning (~ 5 MeV per accreted baryon; 
Schatz et al. 1999) is much larger than that flowing out from deeper in the crust. As a result, the 
temperature at the base of the hydrogen/helium burning shell is determined by the luminosity 
there (for a review, see Bildsten 1998) and is ~ 5 x 10^ K for M ~ 1^~^Mq yr~^ (Brown & Bildsten 
1998; Schatz et al. 1999). The conductivity is much greater in the inner crust (the conductivity 
in the crust increases with density), while the thickness of the inner crust is only a factor of two 
greater than the thickness of the outer crust. For a similar thermal gradient, the inner crust can 
carry a much larger flux. If the change in temperature between the hydrogen/helium burning shell 
and neutron drip is of the same order as the change in the inner crust, than most of the nuclear 
luminosity generated in the crust will flow into the core (Brown &: Bildsten 1998). 



1.2. Outline of this paper 

The remainder of this paper paints in the crude picture just sketched. Sections 2 and 3 develop 
the details of the calculation, which proceeds in two steps. First, the hydrostatic structure of the 
neutron star is computed for different equations of state (§ 2.1) and neutron star masses. These 
hydrostatic models serve as background for the thermal computations, which are discussed in § 3. 
Section 3 also describes the relevant microphysics: the heating by crust reactions and the cooling 
by the crust and core neutrino emissivity (§ 3.1); the reduction of the core neutrino emissivity 
by neutron and proton superfluidity (§ 3.2); and the conductivity of the crust and core (§ 3.3), 
including the effect of impurities. 

Section 4 contains the results of these calculations, which are split into three parts. First, there 
is discussion on the nature of the thermal profile and its dependence on the crust composition and 
the core neutrino emissivity (§ 4.1). Section 4.2 presents some analytical formulae for the crust 
temperature; these formula are used (§ 4.3) to show how the high-accretion rate solutions discussed 
here connect with those at lower accretion rates (e.g., Miralda-Escude et al. 1990; Zdunik et al. 
1992). The melting and refreezing of the crystalline lattice in the inner crust are discussed in § 4.4. 



2. Hydrostatic structure 

Throughout the crust and core, the pressure is supplied by degenerate particles with Fermi 
energies ^ k-QT, and so the equation of state (EOS) scarcely depends on temperature. The crust 
reactions heat the core on a timescale 

( M\CT MCT M 

-;-„ r^ K, ^ (1) 

\muj Ln M En M 
In this equation C is the specific heat per baryon, M/nin is approximately the number of baryons 
in the star, and M/M is the timescale for the mass to increase from accretion. If the heat is stored 



in the electrons (as would be the case if the neutrons were superfluid), the heat content per baryon 
(e.g., Landau & Lifshitz 1980) is CT k, T['^k-BT{k-BT/EY^c) "^ Ejy, where -Ep.e is the electron Fermi 
energy. The heat content per baryon is similar if the neutrons are normal (Lamb & Lamb 1978). 
Equation (1) shows that over timescales necessary to establish a thermal steady-state, the mass of 
the star changes only slightly, and so the hydrostatic equations need not be solved simultaneously 
with the thermal equations. Using a fixed hydrostatic structure simplifies the thermal calculation. 

To calculate the temperature as a function of radius requires integrating the heat transport 
equations over the star. The strong gravitational field modifies the heat flow. In an isothermal star, 
where there is no heat flow, the redshifted temperature is constant, while the proper (as measured 
by a local thermometer) temperature increases as one moves toward the stellar center. Because 
the neutrino emissivity is a strong function of temperature, the thermal transport equations must 
account for gravitational effects. The appropriate equations, solved for stellar mass and EOS, are 
the post-Newtonian stellar structure equations (Thorne 1977) for the radius, gravitational mass, 
potential, and pressure: 



dm p / 2Gm^ ' 

da n \ re 



^(1-^) (3) 
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da Anr'^n \ mc^ J \ rc^ J 



dp Gm p f P \ ( 47rr |)\ / 2Gm,\ 

da Airr'^ n \ pc^ J \ mc^ J \ rc^ J 



(5) 



In these equations the Lagrangian variable a is the total number of baryons inside a sphere of area 
47rr^, and p is the mass density. The potential <I> appears in the time-time component of the metric 
as e^''^ (it governs the redshift of photons and neutrinos; Misner, Thorne, Sz Wheeler 1973) and 
satisfies the boundary condition that at the stellar surface e ''^ \r=R = 1 — 2GA4./ R(? , where Ad 
and 47ri?^ are the total gravitational mass and surface area of the neutron star. 



2.1. Equation of state 

For purposes of calculating the crust EOS, the ashes of hydrogen and helium burning are 
presumed to be pure iron (but see the discussion in § 3.3 on how the composition affects the 
energy release and heat transport). As a mass element is compressed to greater densities and 
pressures, the rising electron Fermi energy triggers a series of electron captures, neutron emissions, 
and pycnonuclear reactions (Sato 1979; Blaes et al. 1990; Haensel & Zdunik 1990b). At any given 
density, only one species is assumed present (see Figure 1) according to the composition calculated 
by Haensel & Zdunik (1990a,b). 



In the outer crust, relativistic degenerate electrons of density Ug = Yen supply the pressure. 
The electron chemical potential is basically the Fermi energy Ep^c = TUeC^ll + (37r^ne)^'^Ag]"^'^, 
where Ag = 386.2 fm is the electron Compton wavelength. I calculate the electron pressure, which 
is approximately ne-E'F,e/4, from the interpolation formula of Paczyhski (1983). The lattice pressure 
is calculated from the ionic free energy, which is a function of 

-irriN , (6) 



kBT V 3 

where n^ is the density of nuclei. I use the fits of Farouki & Hamaguchi (1993) to Monte-Carlo 
simulations of the free energy. (In the crust, the free energy per nucleus is to lowest order just 
the Madelung energy, ~ —O-dTk-sT.) These fits are valid for F > 1, which is always the case for 
the density-temperature regime of interest. Following Farouki &: Hamaguchi (1993), I presume 
the nuclei are crystalline for F > 173. The binding energy of the nuclei are computed from a 
compressible liquid-drop model (Mackie & Baym 1977). This formula accounts for an external 
neutron gas and is therefore applicable at densities greater than neutron drip. The energy density 
and pressure of the neutron gas (which differ from that of an ideal degenerate gas) are also computed 
with this model in the limit of a vanishing proton fraction. 

Summing the pressure contributions from electrons, ions, and neutrons gives the crust EOS, 
which agrees with that of Haensel & Zdunik (1990a) to the accuracy of their table. Figure 2 displays 
this relation, p{n), throughout the inner crust. For reference, the p ex n^'^ relation appropriate for 
an EOS dominated by relativistic, degenerate electrons is also shown {dotted line). Free neutrons 
are present [heavy solid line) for n > 3.6 x 10~^fm~ . As noted by Haensel Sz Zdunik (1990a), for 
n > 0.04 fm~ , the free neutrons provide most of the pressure, and the ionic composition becomes 
less and less important to the EOS. In this regime I use the p{n) fit of Negele & Vautherin (1973). 

At n > 0.1 fm~, the nuclei dissolve into uniform nuclear matter (Pethick et al. 1995). I 
select two sample core equations of state for comparison. The first EOS is a fit (Lai 1994) to 
the AV14+UVII interaction, which is the Argonne V14 potential, with a three-nucleon interaction 
prescribed by the Urbana VII potential (Wiringa, Fiks, &: Fabrocini 1988). The second EOS, 
called AV18+(5v+UIX*, is a Skyrme-type Hamiltonian fit (Akmal, Pandharipande, & Ravenhall 
1998, appendix A) to the Argonne V18 potential, with relativistic boost corrections and the three- 
nucleon interaction UIX* (Akmal et al. 1998). The components of both interactions are neutrons, 
protons, electrons, and, where Ep^^ > ni^c? = 105.66 MeV, muons. I do not consider, for simplicity, 
other possible components (e.g., hyperons or quark matter) in the EOS. To construct a table 
suitable for interpolating n(p), I calculate for each n the proton fraction Yp = n^jn and electron 
fraction Y^ = n^jn from the equations for /3-equilibrium, ^n — p-p = p-e = f^f^i and charge neutrality, 
np = Ue + n^^. Given (Yp, Ye), I then compute the mass density p and pressure p = (P{—p+ndp/dn). 

There have been many attempts to calculate the density of the phase transition from the inner 
crust to the core (see Pethick &: Ravenhall 1995, and references therein). I adopt the following 
approach. The density and pressure of the AV18+(5v+UIX* EOS equal those of Negele & Vautherin 
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Fig. 2. — Pressure as a function of baryon density n in the crust. Each electron capture is visible 
as a discontinuity in n. Where free neutrons are present {heavy solid line) the pressure deviates 
from the piecewise scaling for degenerate, relativistic electrons {dotted line). 



(1973) at n = 0.078 fm^, p = 0.39MeVfni^ . I therefore take this density as the transition from 
crust to core; there is no density discontinuity in this case. For the AV14+UVII EOS, the energy 
density is always greater than that of Negele & Vautherin (1973), and so I choose the maximum 
crust density to be 0.1 fm~' {p = O.GOMeVfm^ ). In this case, there is a substantial density jump 
(from n = 0.1 fm~ to n = 0.13 fm~ ) between crust and core. The choice n = 0.1 fm~ as the 
upper limit for the crust density reflects recent detailed calculations (Pethick et al. 1995) of the 
phase transition. 

For equilibrium crust compositions, it becomes energetically favorable for nuclei to turn inside- 
out in the inner crust and form a phase with bubbles of neutron gas encased in bound nuclear 
matter (Lorenz, Ravenhall, & Pethick 1993; Oyamatsu 1993). Because the charge of nuclei in an 
accreted crust is less than that of the equilibrium composition, it is possible that the nuclei do not 
turn inside-out. The condition for this inversion (see, e.g., Pethick & Ravenhall 1995) is that the 
nuclear radius be more than one-half the Wigner-Seitz radius, (47rn7v/3)~^'^. For the core EOS 
AY18+6v+\]lX*, this ratio at the bottom of the crust is (for Z = 20, A = 100, y„ = 0.8, where 
Yn = nn/n is the neutron fraction) 0.50. For the core EOS AV14-I-UVII, the ratio is 0.54. In the 
absence of a more detailed calculation of the composition, I am unable to say if an accreted crust 
contains non-spherical nuclei and have not explored this possibility. 
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2.2. Construction of models 

With n{p) specified by interpolation from a table, 1 integrate equations (2)-(5) with a fourth- 
order Runge-Kutta integration algorithm (Press et al. 1992). Borrowing a technique used by van 
Riper (1991), the code restricts the stepsize Aa to be always less than some fraction / of the 
reciprocal sum of the radial and baric scale heights. 

Starting from a fixed central pressure, 1 expand the hydrostatic equations about the center a = m = 
r = and integrate outwards until the pressure is less than 3.4 x 10~^ MeV fm~ (corresponding to 
a density n ~ 2.4 x 10^^ fm^ , about a factor of 10 greater than where the helium burning ends; 
Schatz et al. 1999). At this pressure, the radius and mass are constant to within W~^R and 10~^M, 
respectively. The algorithm iteratively adjusts the central pressure until a target gravitational mass 
A4 is reached. The integration steps are then stored for later use in solving the thermal equations. 
For each of the two equations of state, AV14-I-UV11 and AV18-|-(5v-|-UlX*, 1 compute two masses, 
Ai = 1.4Mq and A4 = 1.8Mq; a summary of these four structures is provided in Table 1. 

Despite continuous accretion, the pressure is a good Eulerian coordinate throughout the crust 
(Bildsten et al. 1992; Brown &: Bildsten 1998), and I shall plot the temperature and luminosity 
against it. The mass contained in the crust above a given isobar is to lowest order (from expanding 
eq. [5]; Lorenz et al. 1993) 

AM = m„Aa « -AnR^ 
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where g = GA4{1 + z)/R^ is the gravitational acceleration and (1 + z) = (1 — 2GAA/Rc^)^^''^ is 
the surface redshift. Table 1 lists AM, g, and 1 + z for the two masses and two EOSs considered 
in this paper. 



3. Thermal Structure 

With a hydrostatic structure specified, the luminosity L and temperature T are found by 
solving the entropy and flux equations (Thorne 1977), 

e-2$/c2|-(Le2*/^')-47rr2n(e^-e,)(l-2Gm/rc2)-V2 = q (9) 

e-*/-'K|-('re*/^'K-^(l-2Gm/rc2)-V2 = q. (10) 

or \ ) Airr^ 

Here e^r and e^ are the nuclear heating and neutrino emissivity per baryon, and K is the thermal 
conductivity. I neglect in equation (9) terms arising from compressional heating, as they are of 
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order TAs{M/M) (Fujimoto & Sugimoto 1982), s being the specific entropy, and are negligible 
throughout the degenerate crust and core (Brown &: Bildsten 1998). The physics of the problem 
is contained in the heating ejsf, neutrino cooling e^, and the conductivity K; the following sections 
discuss each in turn. 



3.1. Nuclear heating and neutrino cooling 

As mentioned in § 2.1, the crust electron captures reduce the charge of the nuclei enough 
to trigger pycnonuclear reactions. The rate of pycnonuclear reactions is governed by the supply 
of low-Z nuclei, which is in turn determined by the rate of the preceding electron capture. As a 
result, even though the pycnonuclear reactions are better described at typical crust temperatures as 
strongly screened fusion reactions (Salpeter &; van Horn 1969), they are insensitive to temperature 
and hence not susceptible to a thermal instability. 

To find the overall thermal profile, I do not need to resolve the individual capture layers. 
Rather, I distribute the reaction heat per baryon, E]\f = 1 MeV, over a pressure interval Ap between 
8.7 X 10~*^MeVfm~ and 3.4 x 10~^MeVfm^ (this covers the region where the pycnonuclear 
reactions occur). The total nuclear luminosity is Lj\f = Me^^Ej\f/mn ~ Lyi/200, where La = 
Mc^z/{\ + z) ~ GMM/R is the accretion luminosity (see § 3.4) and Me^* is the accretion rate 
as measured in the crust. The heating term in equation (9) is therefore 

.„.V„.gjMe-/-)(i^|). (U, 

where dp/dr is the Jacobean. 

The cooling terms in equation (9) are evaluated from various fits to microscopic calculations. 
Throughout much of the crust, the dominant neutrino emissivity is neutrino pair bremsstrahlung, 
e~ + {A, Z) -> e- + {A, Z) + uv (Maxwell 1979; Haensel et al. 1996; Itoh et al. 1996). Where the 
ions are crystallized, the bremsstrahlung rates are exponentially suppressed because the separation 
between electron energy bands is of order IMeV ^ k-^T (Pethick & Thorsson 1994, 1997). I 
use the fits of Haensel et al. (1996) for the emissivity where the ions are liquefied and the fits of 
Yakovlev & Kaminker (1996), which include this suppression, where the ions are crystallized. If the 
crust is very impure (see the discussion in § 3.3), then the crust bremsstrahlung will be dominated 
by electron-impurity scattering (Pethick & Thorsson 1997). For T > 10^ K and densities where 
^pe/k-BT > 1, hjJpe ~ 0.056(-E'F,e/lMeV) MeV being the electron plasma frequency, the plasma 
neutrino process (Schinder et al. 1987; Itoh et al. 1996) becomes important. 

This paper assumes a standard core neutrino emissivity, for which modified Urea processes 
(Friman & Maxwell 1979; Yakovlev & Levenfish 1995) dominate. The phase space available for 
scattering is strongly restricted if the nucleons are superfluid and reduces the emissivity roughly 
as exp(— T/Tc) (Yakovlev & Levenfish 1995) where Tc is the superfluid transition temperature. In 
general, both the neutrons and protons must be superfluid to substantially reduce the modified 
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Urea neutrino luminosity. The proton modified Urea braneh {p + p^p + n + e~^ + Vf.) is nearly as 
efficient as the neutron branch (Yakovlev &; Levenfish 1995), and so a slight increase in temperature 
is sufficient to compensate for the suppression of just one of the modified Urea branches. 



3.2. The superfluid transition temperatures 

In the core, both the protons and neutrons are expected to be superfluid over some range of 
densities (Baym, Pethick, & Pines 1969; Hoffberg et al. 1970; Takatsuka & Tamagaki 1993; Amund- 
sen (fe 0stgaard 1985b; Amundsen & 0stgaard 1985a; Elgar0y et al. 1996). At lower densities, the 
neutrons pair in a singlet (^^o) state, but at higher densities the repulsive core of the interaction 
forces the neutrons to pair in a triplet (^^2) state. The protons in the core are expected to be in a 
^^o state. There is at present little agreement on the range of densities for which the protons and 
neutrons are superfluid and on their transition temperatures Tc (for a review, see Pethick et al. 
1995). The early calculation of Hoffberg et al. (1970) for pure neutron matter found peak gap en- 
ergies of 1.6 MeV (singlet) and ^ 5MeV (triplet). They also found that the transition temperature 
remained high over a large range of densities, implying that the entire core would be superfluid. 
More recent calculations of the neutron triplet pairing (Takatsuka & Tamagaki 1993; Amundsen 
& 0stgaard 1985b) find a lower transition temperature k-^T^^^ ^0.1 MeV. In addition, the range 
of densities for which pairing occurs is restricted. Elgar0y et al. (1996), using a meson-exchange 
model for /3-stable matter, found that maximum gap energy was ~ 0.018 MeV and that the range 
of densities was restricted to n < 0.13 fm~, so that almost all of the core would be normal. In 
all cases the proton pairing gaps are somewhat larger, with k-^T^^^ ^ 1.0 MeV, and extend to 
densities several times the saturation density, Ug = 0.16 fm~. 

Although none of the published microscopic calculations of the critical temperature Tc has 
presented a convenient fitting formula in terms of density, the critical temperatures are roughly 
quadratic functions of the Fermi wavevector ks^p} = (3vr^nr„p|)^'^. I therefore use 



Tc{k) = T, 



cO 



{k-h^^ 



^0) 



(12) 



(Afc/2)2 

as the functional form of Tc for the proton ^5*0, neutron ^5o, and neutron ^^2 states. The parameters 
Tco, A^O) ^-iid Afc are chosen (see Table 2) to approximate the transition temperature of Amundsen & 
0stgaard (1985b) for the neutron ^P2 state and the transition temperature of Amundsen &; 0stgaard 
(1985a) for the proton and neutron singlet states. This choice of Tc reflects the calculations of 
Takatsuka & Tamagaki (1993) as well. 

Figure 3 displays Tc for the proton ^Sq {solid lines) and neutron ^^2 {dashed lines) pairing as a 
function of density, for both equations of state. While the neutron critical temperatures are roughly 
identical, the proton critical temperature for AV18+5v+UIX* vanishes at lower n. This cutoff is 
because AV18+(5v+UIX* has a higher proton fraction, at a given n, than AV14+UVII. For each 
EOS, I show {arrows) the central densities of neutron stars of gravitational masses A4 = IAMq 
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Table 1. Hydrostatic structures 



model 


M 


M 


R 


l + z 




g 


Pcrust 


AM 




(Mo) 


(Mo) 


(km) 




(10^4 


cms^^) 


(MeVfm^^) 


(Mo) 


Ml.4-14^ 


1.58 


1.4 


10.7 


1.28 


2.08 




0.60 


0.025 


Ml. 8-14^ 


2.12 


1.8 


10.4 


1.43 


3.13 






0.017 


Ml.4-18*' 


1.56 


1.4 


11.6 


1.25 


1.74 




0.39 


0.021 


Ml.S-lS*' 


2.09 


1.8 


11.3 


1.38 


2.61 






0.015 



''EOS is a fit (Lai 1994) to AV14+UVII (Wiringa et al. 1988) 
i^EOS is AV18+(5v+UIX* (Akmal et al. 1998, appendix A) 

Note. — M is the number of baryons in the star, multiplied by m„ = 1.67 x 10~^^ g; 
M. is the gravitational mass of the star; Pcrust is the pressure at the base of the crust; 
and AM is the number of baryons in the crust multiplied by nin- 



Table 2. Parameters used to calculate superfluid transition temperature (eq. [12]). 



State Tc ko Ak 



MeV) (fm^^) (fm" 



proton ^So 0.345 0.7 1.0 

neutron ^Sq 0.802 0.7 1.2 

neutron ^P2 0.0776 2.0 1.6 
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Fig. 3. — Superfluid transition temperature T^ (see eq. [12]) as a function of density, for both the 
proton ^5*0 {solid lines) and neutron ^P2 {dashed lines) pairing. Each panel is for a different core 
EOS, as indicated in the legend. The arrows mark the central densities, for each EOS, of a neutron 
star of gravitational mass A4 = IAMq and Ai = 1.8Mq. 

and Ai = I.SMq. At core temperatures ~ 5 x 10^ K, each of the four EOS/mass combinations 
has normal protons, neutrons, or both in some part of the core. As a result, the modified Urea 
processes still play an important role in the neutron star's thermal balance^. 



3.3. Heat transport 

Throughout the crust, relativistic electrons transport the heat. In the relaxation-time approx- 
imation, the conductivity is (e.g., Ziman 1972), 



K 



■K k^Tn^ 



-fee 



mz 



(13) 



Here m* = E-p^^/f? is the effective electron mass, and r is the relaxation time. Where the ions 



are crystallized, l/r 
In these formulae, r^, 



1/ree + I/tcQ + 1/rep, and where they are hquefied, l/r = l/xee + 1/rej. 

TeQ, Tgp, and Tei are respectively the relaxation times for electron-electron 
(Urpin &: Yakovlev 1980; Potekhin, Chabrier, & Yakovlev 1997), electron-impurity (Itoh & Ko- 
hyama 1993), electron-phonon (Baiko & Yakovlev 1995), and electron- ion (Yakovlev & Urpin 1980) 
scattering. Electron-electron scattering is typically negligible over much of the crust because the 
strong degeneracy of the electrons restricts the available phase space. 



Recently, neutrino emission from the formation and destruction of Cooper pairs has received renewed interest 
(see Yakovlev, Kaminker, & Levenfish 1999, and references therein). For temperatures near the superfluid transition 
temperature, the neutrino emissivity is enhanced over that of modified Urea processes. I have not included this 
emissivity here; for an impure crust, this omission is not critical (see § 4.2). 
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The electron-impurity relaxation time is (Yakovlev & Urpin 1980) 

Here pp and vp are the momentum and velocity of an electron at the Fermi surface and AgQ ~ 2 
(Yakovlev & Urpin 1980) is the logarithmic Coulomb factor. The concentration of impurities enters 
through Q, which for a large number of ion species is (Itoh &: Kohyama 1993) 

Q=—Tn,{Z,-{Z)f. (15) 

nN ^-^ 

Here n^ and Zi are the density and charge number of the ith species, nj^ = Y2^ ni is the total ionic 
density, and {Z) = n^ ^^ njZj is the mean ionic charge number. 

There are several contributions to the impurities in the crust. First, within each electron 
capture layer there are at least two species present. The electron capture from an even-even 
to an odd-odd nucleus is immediately followed by a second electron capture to a lower energy 
even-even nucleus (Haensel & Zdunik 1990b). Within each capture layer the impurity parameter 
is then Q = Anznz-2/{fiz + nz~2f' l£ !> with nz and nz-2 denoting the densities of the two 
species. Although Haensel &: Zdunik (1990a,b) treated the electron captures as sharp transitions 
in the crust, in actuality the layers have a finite thickness set by competition between the flow 
timescale and the (weak) electron capture timescale (Bisnovatyi-Kogan & Chechetkin 1979; Blaes 
et al. 1990; Bildsten & Gumming 1998). The zero-temperature electron capture rate is proportional 
to {Ep^c — 6)^, where 6 is the reaction threshold. Where degenerate electrons supply the pressure, 
£'f,c must increase with depth. The layers are thin in Ep^^, and because p oc Ep ^, the layers are 
geometrically thin as well. For densities greater than neutron drip, however, the electrons no longer 
support the crust, and -Ep.e need not increase with depth. In fact, £'F,e is actually less following 
an electron capture layer if the interface between layers of different composition is treated as a 
infinitely thin plane. The layers, although thin with respect to Ep^e, are then geometrically thick. 
In actuality, thermal broadening of the electron Fermi surface causes many of the captures to occur 
pre-threshold, and the capture layers are thickened to nearly the width between layers (Ushomirsky 
et al. 1999). Should the capture layers overlap, then Q in the mixed layer can become larger than 
unity. 

The impurities manufactured within the capture layers are probably a small perturbation 
compared to those already present in the mixture entering the top of the crust. Schatz et al. 
(1999) found that Q ~ 100 immediately following the end of stable hydrogen burning. An accurate 
assessment of Q, throughout the crust, requires evolving the composition of an accreted fluid 
element on its journey through the crust. This task is beyond the scope of this initial survey, and 
I instead set upper and lower bounds on the conductivity. The upper bound to the conductivity 
is that of a pure crystal (electron-phonon scattering). To set the lower bound, first note that 
electron-impurity scattering dominates the conductivity wherever TeQ < T(.p, with Tgp being the 
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electron-phonon relaxation time (Baiko &; Yakovlev 1995) 

Here a is the fine structure constant, and Agp ~ 13 comes from integrating over the phonon 
spectrum. Equations (14) and (16) imply that for Q > 0.66(30 MeV/^F,c)(^/26)(A;BT/0.05MeV) 
electron-impurity scattering determines the thermal conductivity in the crust. If the reactions in the 
crust do not significantly reduce Q from its large value at the base of the hydrogen/helium burning 
shell, then the heat transport in the crust is controlled by electron-impurity scattering rather than 
electron-phonon scattering. 

If Q is very large ('^ Z^), then the impurity relaxation time is roughly that of electron- ion 
scattering for a pure crystal (Yakovlev & Urpin 1980), 

- - "^"^ A-', (17) 



with Aei = \n[{2TxZ/2>y-'^ ^Jl^b+'^]T] — 1. Basically, the phonon spectrum is extremely disordered 
in this case. I therefore set a lower limit to the conductivity by using electron-ion scattering, i.e., 
by treating the ions as if they were liquefied. For consistency, I also use the liquid-state neutrino 
bremsstrahlung emissivity (Haensel et al. 1996) in conjunction with the electron-ion conductivity. 
The impurities in the crust reduce the conductivity but increase the neutrino emissivity. 

In the core, heat is mostly carried by electrons, with neutrons contributing if they are normal 
(Flowers & Itoh 1979). I neglect here the neutron conductivity. This is a good approximation, as 
the core is practically isothermal (see § 4). In evaluating the electron-proton scattering terms, I 
used an effective proton mass m* = 0.7mp. The proton superfluidity both reduces the screening 
(increases the scattering) and reduces the proton scattering phase space (suppresses the scattering) ; 
I take these factors into account using the fits of Gnedin & Yakovlev (1995). 



3.4. Boundary conditions and method of solution 

The first boundary condition is simply L\r=o = 0. For a fully self-consistent solution, the correct 
second boundary condition is a relation L(T)\r=R, usually obtained from a separate photospheric 
calculation. This is unnecessary, however, when the hydrogen and helium burn steadily. The large 
energy release from this burning determines the temperature in the outer atmosphere, so that the 
temperature at the base of the hydrogen/helium burning shell is a function only of M, M, and R 
and may be calculated independently. In § 4, I show that the luminosity flowing outwards from the 
crust is in fact much smaller than that generated by the hydrogen/helium burning. 

The second boundary condition, then, is T\r=R = To{m), where I take To from Schatz et al. 
(1999). Here rh is the accretion rate per unit area; the fiducial rate used by Schatz et al. (1999) is 
the Eddington rate appropriate for a Newtonian star of A4 = 1.4M0 and R = 10 km accreting a 
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solar composition plasma. Numerically, this rate is rriE = 8.8 x lO^gcm^^s^"^ and is an excellent 
approximation to the lowest local accretion rate at which the hydrogen/helium burning is stable 
(see, e.g., Bildsten 1998). At m = ifiE, the temperature at the base of the hydrogen/helium burning 
shell is To = 5 X 10*^ K. Although the Newtonian surface gravity used by Schatz et al. (1999) is less 
than the values used here, To is relatively insensitive to g (To oc g^''^; Bildsten 1998), so I do not 
adjust it for each model. This is not critical, as the thermal profile in the crust is insensitive to To 
(see § 4). 

The accretion rate enters equations (9) and (10) through e^, which is scaled to the accretion 
luminosity La- Because To is only a function of m, I use the same rh for each model; the luminosity 
from this accretion is then different for each EOS and mass. The global accretion rate, as measured 
by an observer infinitely far away, is (Ayasli & Joss 1982) M = 47ri?^m/(l + z), and the luminosity 
is 

La = ^—Mc^ = - — ^-Tr47r i^^mc^. (18) 

l + Z (1 + Z)2 ^ ^ 

For the fiducial local accretion rate m^;, model Ml. 4-14 has La = L°^ = 1.95 x 10^^ ergs~^; Ml. 8-14 
has L\ = 2.28 x lO^^ergs"^; Ml.4-18, 2.12 x lO^^ergs-^; and Ml.8-18, 2.51 x lO^^ergs'^. To 
solve the thermal structure, equations (9) and (10) are finite-differenced onto the mesh defined 
by the integration of equations (2)-(5). An initial guess is constructed by fixing the temperature 
throughout the star to To and integrating equation (9) from L\r=Q = 0. This trial guess is then 
iteratively refined by a relaxation technique (Press et al. 1992). The resolution of the mesh was 
tested by computing models with step fractions (see eq. [7]) / = 0.05 and / = 0.02. 



4. Results 

4.1. The influence of the microphysics on the thermal profile 

As promised in section 3.3, I survey the uncertainties in the thermal conductivity by solving 
the thermal structure (eq. [9] and eq. [10]) with the conductivity alternately set by electron-phonon 
scattering and electron-ion scattering. Figure 4 shows the thermal profiles for these two cases with 
different degrees of superfluidity: strong (both neutrons and protons are superfluid with T^.'^T 
throughout the core; top panel) ^ moderate (corresponding to the parameters in Table 2; middle 
panel), and nonexistent [Tc = OK for both neutrons and protons; bottom panel). The hydrostatic 
structure used in this plot is model Ml.4-18 (see Table 1). I plot the proper temperature (i.e., the 
temperature a local thermometer would measure) because it controls the conductivity and neutrino 
emissivity. 

The lower conductivity and enhanced bremsstrahlung from electron- ion scattering {dotted 
lines) produce a greater temperature variation throughout the crust than if the conductivity and 
bremsstrahlung were determined by electron-phonon scattering (solid lines). In the inner crust, 
the electron-phonon conductivity (Fig. 5, solid line) is an order of magnitude greater than the 
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Fig. 4. — Proper temperature, in units of 10^ K, as a function of pressure. The hydrostatic structure 
is Ml. 4-18 (see Table 1). Two cases are compared, with the conductivity dominated either by 
electron-phonon scattering {solid lines) or by electron-ion scattering {dotted lines). The three 
panels, from top to bottom, show the variation of the thermal profile with core superfluidity. In 
the top panel, the core neutrino emission is completely suppressed: both protons and neutrons 
are superfluid throughout the core. The middle panel shows the effects of moderate superfluidity 
(see Table 2), whereas in the bottom panel there is no superfluidity. Note that for a low thermal 
conductivity, the maximum crust temperature (~ 8 x 10^ K) varies only slightly with the core 
temperature. The pressures of neutron drip and the crust-core boundary (see Table 1) are marked 
along the bottom axis of each plot with arrows. 
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Fig. 5. — Thermal conductivity in the accreted crust. The thermal structure corresponds to the 
cases shown in the middle panel of Figure 4. Electron-phonon is denoted with a solid line, while 
electron-ion is denoted by a dotted line. Note that there is an order of magnitude difference between 
the two cases in the inner crust. 



electron-ion conductivity (Fig. 5, dotted line). As a result, the thermal gradient in a locally pure 
crust is very small, so that the inner crust temperature is not appreciably different from that of 
the core. A striking result for this case is that the peak crust temperature depends only weakly 
on the core temperature. This is a consequence of the requirement that a large thermal gradient is 
needed to carry the flux in the inner crust when the conductivity is reduced. 

The relative amounts of neutrino emission from the crust and core are displayed in Figure 6, 
which shows the luminosity measured by an observer at infinite distance and scaled to L^. The 
lines and panels correspond to those of Figure 4, and for reference the region where nuclear heating 
occurs (see § 3.1) is denoted with boldfaced lines. A negative luminosity indicates that the heat 
flow is inward. When Tc ^ T throughout the core for both neutrons and protons (top panel), all 
of the neutrinos are emitted from the crust, and the luminosity is zero throughout the core. For 
moderate superfluidity (Table 2; middle panel), some neutrino emission occurs in the crust when the 
conductivity and crust bremsstrahlung are dominated by electron-ion scattering {dotted line), but 
the bulk of the neutrinos are emitted in the innermost core, where the protons are normal. If there 
is no superfluidity, then neutrino emission occurs throughout the core (bottom panel). The decrease 
in L at pressures < 10^'^ MeV fm~ marks where the plasma neutrino emissivity dominates. 

The core temperatures of the middle and bottom panels of Figure 4 are similar because the 
modified Urea emissivity is strongly temperature sensitive, so that the core temperature need only 
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Fig. 6. — Luminosity, as measured by an infinitely distant observer and in units of the total crust 
nuclear luminosity Ln, plotted as a function of pressure. Negative values indicate an inwardly 
directed flux. The panels correspond to those in Figure 4, and both high conductivity (electron- 
phonon scattering; solid lines) and low conductivity (electron-ion scattering; dotted lines) cases are 
shown. The region where nuclear heating occurs is indicated with boldfaced lines. 
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Fig. 7. — Proper temperature, in units of 10^ K, as a function of pressure for two different core 
equations of state, AV14+UVII {solid lines) and AV18+5v+UIX* {dotted lines), and two different 
gravitational masses, Ai = 1.4M0 {top panel) and A4 = 1.8Mq {bottom panel). The superfluidity 
is as described in Table 2, and the conductivity is set by election-ion scattering throughout the 
crust. The pressures of the crust-core boundaries for the two EOS is given in Table 1. 



slightly adjust to compensate for a reduced normal core fraction. To demonstrate this further. 
Figures 7 and 8 show, as functions of pressure, the proper temperatures and luminosities for the 
four models in Table 1 with conductivity set by electron-ion scattering. The equations of state 
AV14+UVII {solid lines) and AV18-F(5v-MJIX* {dotted lines) are compared for M = IAMq {top 
panel) and A4 = I.8M0 {bottom panel). For A4 = 1.8Mq the two equations of state have similar 
thermal profiles because both have normal protons and neutrons in at least some fraction of the 
core (see Figure 3). In contrast, for M = 1.4M0, only Ml. 4-18 {top panel, dotted line) has normal 
protons in its innermost core, and so its core temperature is quite cooler than that of Ml. 4-14 {top 
panel, solid line). As a result, Ml. 4-14 has a stronger neutrino emission from the crust (Figure 8, 
top panel, solid line). As in Figure 6, the nuclear heating region is denoted with boldfaced lines. 

A generic feature of these solutions is that almost all of the nuclear heat released in the inner 
crust flows inward and is balanced by neutrino cooling from either the crust or core. Only a small 
amount (< 5%) of Ljv is conducted to the surface. As a result, the temperature in the inner 
crust and core is set by the processes in the inner crust. In particular, the temperature of the 



22 




•e- 
0) 



-1 



-0.5 

-1 



I I I I I M M I I I I I I I 



Jl = 1.8 M, 



neutron drip 



± 



I I I I I 



-8 -6 -4-2 2 

log(p/MeV fm"2) 

Fig. 8. — Luminosities corresponding to the thermal profiles in Fig. 7, with negative values indi- 
cating that the flux is inwardly directed. The higher neutrino emissivity from the modified Urea 
proton branch is evident for model Ml. 4-18 {top panel, dotted line): the AV18-|-5v-|-UlX* EOS has 
a larger proton fraction Yp and the protons in the inner portion of the core are normal (see Fig. 3) . 
Nuclear heating occurs in the region denoted with boldfaced lines. 
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Fig. 9. — Proper temperature, in units of 10^ K, as a function of pressure. A family of solutions for 
different To, the temperature at the top of the crust, is shown, both for a thermal conductivity set by 
electron-ion scattering {top panel) and electron-phonon scattering (bottom panel). The temperature 
sensitivity of modified Urea processes (superfluidity as described in Table 2) keeps the base of the 
crust at a temperature ~ 4.3 x 10^ K in both cases. 

inner crust is nearly independent of the temperature in the hydrogen/helium burning shell. This is 
explicitly shown in Figure 9, where I plot the thermal structure for model Ml. 4-18 but with the 
outer boundary temperature allowed to vary. The luminosity is fixed at L^ = 2.12 x 10^^ ergs~^. 
The top panel shows the case of electron-ion conductivity; the bottom, electron-phonon. 



4.2. Analytical expressions for the crust temperature 

The numerical results presented above can be easily understood by considering a crude analyt- 
ical calculation of the thermal structure. The approach is similar to that of Hernquist & Applegate 
(1984), with three exceptions: I include heating from the crust reactions, I fit the pressure as a func- 
tion of density, rather than presume a degenerate relativistic EOS, and I assume the conductivity 
is determined by electron-ion, rather than electron-phonon scattering. 



In the crust, the surface gravity and radius are roughly constant, and the differential expressions 
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for the radius (eq. [2]), pressure (eq. [5]) and flux (eq. [10]) can be combined into the plane-parahel 

Newtonian equation 

dT L 

To construct this analytical model, I consider the heating and cooling emissivities to be (5-functions. 
Since the flux is then constant between points where these sources or sinks reside, I may integrate 
equation (19) piecewise between these points, with L stepping discontinuously at each point (Brown 
& Bildsten 1998). 

Integrating equation (19) requires a relation p{p). I approximate the mass density p by rriun 
and flt the pressure with power-laws in both the electron-dominated and neutron-dominated regions. 



P = 



2.67 X 103Vl2^^dynecm-^ pi2 < 0.66 

A n^ .. ir.29 AA2 J -2 _ ^ o r. 



4.97 



J^u />i2 uynecm , pi2 ^ u.uu, , . 

10^^12''^ dyne cm-^ P12 > 8.9. ^ ' 

Here pi2 = p/lO^^ gcm^^, and the error in p is less than 9% and 5% for the two density regimes 
respectively^. The exponent in the electron-dominated regime is less than 4/3 because the fit 
accounts for the decrease in Yg by electron captures. For densities just above neutron drip (0.6607 < 
P12 < 8.913), the pressure cannot be fit by a simple power-law in density (see Figure 2). This region 
is where most of the heat is released, and so for simplicity I presume the region to be isothermal 
and place the heating (5- function inside it. 

Inserting the expression for the electron- ion scattering frequency (eq. [14]) into the expression 
for the thermal conductivity (eq. [13]) and expanding, I have 

K « -( — 1 a-^h''^-^[^\ klT 




pJ^Tgergcm^^s^^K^^ (21) 

In this expression, I set the Coulomb logarithm Agi to unity and use the shorthand Tg = T/IO^ K. 
The composition enters through Ye /Z] for densities less than neutron drip I use Ye jZ = 

0.0298, appropriate for a pure iron composition, and for densities greater than neutron drip I use 

1/3 
Ye /Z = 0.0167, as follows from the last entry of Table 2 in Haensel & Zdunik (1990b). (For 

impurity scattering, Ye /Z ^ ZYe /Q.) Using equations (20) and (21), I integrate equation (19) 

from the top of the crust to neutron drip, 10~^ < pi2 < 0.66, and over the inner crust to the core, 

8.913 < pi2 < 166, to obtain 

r|(pi2) = T^D^g - 1.30Lo,35^ro'5r4'3 (pr2°-°'° - 1-03) (22) 

for P12 < 0.66 and 

Tiipu) = T^D,9 + 0.33Li^35RTo9il3 {pu'' " l-2l) (23) 



'in this section I use cgs units for easy comparison with the astrophysical literature. 
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for pi2 > 8.9. In these equations Tnd is the temperature at neutron drip (presumed constant for 
0.66 < pi2 < 8.9), Rio = i?/10km, (714.3 = g/10^'^-^ cms~^, and ^0,35 and -Lj,35 are the luminosities, 
in units of 10^^ erg s^^, for pi2 < 0.66 and pi2 > 8.9, respectively. Both Lq and L, are signed: 
they are positive if the flux is directed outwards and negative if directed inwards. Notice that the 
coefficient of Lo,35 in equation (22) is an order of magnitude larger than the coefficient of Lj^35 in 
equation (23). This disparity reflects that the inner crust requires a much smaller thermal gradient 
than the outer crust to carry a given flux. 

To solve for the thermal structure, I also require that the luminosity flowing away from the 
crust heat source is 

Ln = Lo- U, (24) 

and that the core neutrino luminosity balances the heat conducted into the core, 

Li + L^(rcorc) = 0. (25) 

Evaluating equation (22) at pi2 = 10^^ and equation (23) at pi2 = 166, and using equations (24) 
and (25) to replace Lq and Li with L^ and Ly, I obtain an equation for the core temperature, 

^core.g = ^0,9 + Rio9iL [l-26L;v,35 - 1.38L,,35 (Tcorcg)] ■ (26) 

Here To = T\p^^^iQ-5. Solving equation (26) for a modified Urea luminosity Ljy^35 k, b x IQ'^T^^^^^ 
and Ln ~ 1.07 x 10^*^ ergs"^ (m/rhs) gives Tcorc,9 = 0.34 and Li,,35 = 0.92LAr,35, i.e., 92% of the 
heat generated in the crust flows into the core. This compares reasonably well with the numerical 
calculation without core superfluidity (Figures 4 and 6, bottom panels). In that case, the proper 
temperature at the crust bottom is 3.1 x 10^ K, the luminosity flowing out the top of the crust 
is 0.04LAr, and the luminosity flowing into the core is O.SL^. Neutrino emission from the crust 
balances the remainder of L]\f. 

Substituting Tcore,9 and L^^^d into equation (23), I find that TnDjQ = 1-1, which is an overes- 
timation of the maximum crust temperature. This is a consequence of the "two-zone" treatment, 
which puts all of the neutrino cooling in the core. Still, the qualitative features of the numerical 
solutions are reproduced. In equation (23), the increase in temperature from core to neutron drip 
(second term, right-hand side) is much larger than T^^^^g- As a result, changing Tcore,9 has only 
a small effect on Tnd,9- Even if the direct Urea process were to operate and cool the core to 
Tcorejg *^ 1) the temperature around neutron drip will remain high. At high accretion rates, the 
temperature around neutron drip, for a very impure crust, is primarily determined by the ability of 
the inner crust to carry the nuclear luminosity inward and not so much by the efficiency of core 
neutrino cooling. 



4.3. Accretion at higher and lower rates 

As the crust temperature increases, crust neutrino bremsstrahlung and the plasma neutrino 
process become increasingly important. At the higher accretion rate, the brighter crust neutrino 
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Fig. 10. — Proper temperature, in units of 10^ K {top panel), and luminosity, scaled to the nuclear 
luminosity, as measured by an observer at infinity (bottom panel) for accretion luminosities L°^ = 
2.12 X 10^^ ergs~^ {solid lines) and 5L^ = 1.06 x 10^^ ergs""^ {dotted lines). At the higher accretion 
rate, the cooling from crust bremsstrahlung and plasma neutrino emission balance the heating 
{boldfaced lines) on the spot, so that the luminosity flowing out the top of the crust is nearly zero. 



luminosity balances the nuclear heating "on the spot." Figure 10 compares proper temperature {top 
panel) and scaled luminosity {bottom panel), as measured by an observer at infinite distance, for 
model Ml. 4-18 accreting at rriE {La = L\ = 2.12 x 10^^ ergs~^; solid line) and fyiriE {La = ^L\ = 
1.06 X lO'^^ergs"^; dotted line). The conductivity in both cases is set by electron-ion scattering. 
As the crust neutrino cooling increases, a smaller fraction of Ljv flows outward from the top of the 
crust. 

At lower accretion rates, the change in temperature over the inner crust becomes smaller 
relative to the core temperature. The crust becomes more nearly isothermal and hence more 
sensitive to the temperatures at its boundaries (cf. Miralda-Escude et al. 1990; Zdunik et al. 1992). 
From equation (23), the temperature increase over the inner crust is < 0.5Tcorc for La < 0.06L^, 
assuming that Lj = L^. To demonstrate this. Figure 11 displays, for an accretion luminosity La = 
O.OIL^ = 2.12 X 10^^ ergs~^, the proper temperature and luminosity. The hydrostatic structure is 
Ml. 4-18, the same as in Figure 9. The top panel is for a conductivity set by electron-ion scattering; 
the bottom panel is for a conductivity set by electron-phonon scattering. Solutions for several To 



27 




1^ 



^ 0.25 - 



n drip crust/core 

I I I I ^ I I I I ^ I I 



-6 -4 -2 

log(p/MeV fm"3) 



Fig. 11. — Temperature, in units of 10^ K, as a function of pressure for an accretion luminosity 1/100 
of the fiducial rate plotted in Fig. 9. Both the high-conductivity (electron-phonon scattering; bottom 
panels) and low-conductivity (electron-ion scattering; top panels) cases are considered. Unlike 
Fig. 9, the temperature in the inner crust and core depends more on the temperature at the top of 
the crust. The sharp kink in the thermal profiles {bottom panel) is where the ions crystallize. 



are shown; the range of values are reduced from those used in Figure 9 by (m/m^;)^''^ = 0.01^''^, 
which is roughly how the temperature at the base of a hydrogen/helium burning shell scales with 
accretion rate (Schatz et al. 1999). Of course, the hydrogen and helium ignition is unstable in 
an envelope this cold (see Bildsten 1998, and references therein), and so To is determined by the 
compression of matter in the atmosphere and by the flux flowing out the top of the crust. As To 
is reduced, more and more of the heat generated in the crust flows outwards rather than into the 
core. As found by Zdunik et al. (1992), an enhanced core neutrino emissivity will drastically lower 
the crust temperature for low accretion rates. 

To illustrate how the crust temperature changes with the temperature in the hydrogen/helium 
burning region (To), I compute the derivative dTcmst / dTo , where Tcrust is the temperature at the 
centroid of the heat-producing region, p = 0.017 MeVfm^. Figure 12 displays dTcmst/dTo as 
a function of To for five different accretion rates: La/L\ = 0.01 {hollow triangles), 0.03 [filled 
triangles), 0.1 {hollow squares), 0.3 {filled squares), and 1.0 {asterisks). When the conductivity is 
low (electron-ion scattering; top panel), Tcrust is generally less sensitive to To than when electron- 
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Fig. 12. — Derivative of the crust temperature (at the centroid of the heat-producing region, 
p = 0.017 MeV fm^ ) with respect to the temperature at the top of the crust (in units of 10^ K) for 
Lj^/L°^ = 0.01 (hollow triangles), 0.03 (filled triangles), 0.1 (^hollow squares), 0.3 (filled squares), 
and 1.0 (asterisks). The top panel is for a low conductivity (electron-ion scattering) crust, while 
the bottom is for a high conductivity (electron-phonon scattering) crust. 

phonon scattering sets the heat transport (bottom panel). The derivative (at a given accretion rate) 
is largest when To = Tcrust! this peak is evident in the top panel for La/L\ = 0.01 and in the bottom 
panel for La/L\ = 0.03. The rapid rise of dTcmst/dTo in the bottom panel for La/L\ = 0.01 is 
because the neutrino cooling in the crust and core goes to zero, so that all of the heat generated 
in the crust flows outwards (cf. Figure 11, bottom panel). In addition, the crust in the entire 
region considered is also crystalline, which reduces dT/dp. In general, for M > 10^^Mq/yt:~^, the 
temperature in the crust becomes independent of the temperature in the atmosphere and upper 
ocean of the neutron star. 



4.4. Crust melting 

An interesting possibility for a rapidly accreting neutron star is that its crust may melt. This 
happens wherever F < 170, where the exact value is uncertain (for a review, see Ichimaru 1982). 
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Fig. 13. — Ion coupling parameter F as a function of pressure in the neutron star crust, for Ml. 4-18 
(top panel) and Ml. 8-18 {bottom panel), as described in Table 1. The conductivity is alternately 
dominated by electron-phonon scattering (solid lines) and electron- ion scattering (dotted lines). 1 
also display the melting criterion F = 173 (heavy solid lines). 



Since 1 use the formulation of Farouki & Hamaguchi (1993) to calculate the ionic free energy, 1 also 
adopt their melting value, F^ = 173. The crust reactions reduce Z and heat the crust; both of 
these effects decrease F, as shown in Figure 13 for models Ml.4-18 (top panel) and Ml. 8-18 (bottom 
panel), with each model accreting at its fiducial rate. In both cases, the core superfluidity is as 
described in Table 2. For a low thermal conductivity (dotted lines), the crust melts in a series of 
layers. The jaggedness of F is because of the pycnonuclear reactions. Each one doubles Z and 
halves nj\f, so that F increases by 2^'^ and the crust refreezes. Electron captures then decrease Z 
and F until the crust melts again. 

As a consequence of this melting and freezing, the crust resembles a layer cake at densities 
greater than neutron drip. Figure 14 shows the nuclear charge Zm (thin line), below which the 
ions are liquid, along with the Z of the nuclei present (thick lines) according to Haensel Sz Zdunik 
(1990a). The thermal structure is the same as plotted in Figure 13, top panel, dotted line. F 
increases with density (or equivalently, Zm decreases), and so the naive expectation is a sharp 
transition from an ionic ocean to a crust. The Fermi energy also increases with density, however. 
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Fig. 14. — Nuclear charge Zm {thin solid line) such that a pure lattice composed of ions of charge 
Z < Zm melts, as a function of density n. This solution is for the same parameters as Fig. 13, top 
panel, dotted line. The ions are liquid in the region below this curve and crystalline in the region 
above it. Also plotted are the Z of the nuclei present at each depth (heavy lines), according to 
Haensel & Zdunik (1990a). 



and the ensuing decrease in Z from electron captures offsets the rise in F: both Zm and Z decrease 
together. The melting strongly depends on composition: an increment of Z by 2-3 is enough to 
keep the crust crystalline throughout. 

Notice from Figure 13 that there is no crustal melting if the conductivity is solely determined 
by electron-phonon scattering, for models Ml. 4-18 and Ml. 8-18 (solid lines). For model Ml. 4-14, 
the crust melts even if the conductivity is set by electron-phonon scattering. For larger masses 
(models Ml. 8-14 and Ml. 8-18), a high impurity concentration is needed to ensure melting. This is 
a consequence of the stronger core neutrino cooling holding the crust at a slightly lower temperature 
(cf. Figures 7 and 8). Of course, when the impurity concentration is high, the single-species 
calculation of F (eq. [6]) is no longer applicable. Calculations for binary-ionic mixtures (e.g., 
Segretain & Chabrier 1993) show that the melting temperature is lowered below that of the pure 
phases. While the phase diagram of a plasma composed of a large number of species has not 
been calculated, it is likely that the Tm of an impure crust is lower than that assumed here. 
This strengthens the contention that the impure crust of a rapidly accreting neutron star contains 
melted layers, provided that the Z used here (from Haensel &: Zdunik 1990a) is roughly the average 
charge of the nuclei actually present. A self-consistent calculation of the crust composition, and the 
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resulting phase diagram, is required to conclusively determine if layer cake melting actually occurs. 



5. Summary and concluding remarks 

There are three main conclusions presented in this work. First, for neutron stars accreting 
rapidly enough for the accreted hydrogen and helium to burn stably, most of the heat released 
in the crust flows into the core. As a result, the thermal profile in the inner crust is nearly 
independent of the temperature at the top of the crust. Second, if the crust lattice is very impure, 
there is a maximum in temperature at densities greater than neutron drip, where the heating 
occurs. The peak temperature in the crust in this case is set by the ability of the crust to carry the 
generated nuclear luminosity inward from the reaction shell and is relatively insensitive to the core 
temperature. Third, heating the inner crust to temperatures ~ 8 x 10^ K might melt the crust in 
thin layers where electron captures have reduced the ionic charge. 

There are several consequences of these results. Because a fluid layer does not support shear 
stress, the strain in the crust must vanish in these melt layers. This will limit the quadrupole that 
can be induced by thermal perturbations to the electron capture rate (Bildsten 1998) if these cap- 
tures occur above the melt layer. In addition, the fluid layers can dissipate rotational energy, either 
through hydro dynamical or magnetohydrodynamical processes, and thus contribute to balancing 
the accretion torque acting on the stellar surface. The electrical conductivity of an accreted crust 
is reduced, both because of crust heating (Urpin & Geppert 1995; Geppert & Urpin 1994) and 
because of crust impurities (Brown & Bildsten 1998). If the crust is as impure as considered here, 
the timescale for Ohmic decay over a pressure scaleheight is much less (by a factor of 100) than the 
flow timescale, for much of the crust. As a result, the inward advection of magnetic flux (Konar 
& Bhattacharya 1997) is reduced in importance. Thermomagnetic effects, such as current drift 
(Geppert & Urpin 1994) and the battery effect (e.g., Blandford, Applegate, & Hernquist 1983), will 
be comparatively more important, however, because of the greater thermal gradient. 

In recent years, attention has been given to other, more efficient, cooling mechanisms. The 
direct Urea process can operate if the proton fraction is greater than 0.148 (Lattimer et al. 1991) or 
if hyperons are present (Prakash et al. 1992). Other exotic mechanisms may be possible, including 
pion condensates (Umeda et al. 1994), kaon condensates (Brown et al. 1988), or quark matter 
(Iwamoto 1982). The exotic mechanisms have the same temperature dependence as the direct Urea 
(oc T®) but are weaker. Although none of the hydrostatic structures considered in this paper has 
an interior proton fraction large enough to activate the direct Urea, some form of enhanced cooling 
could operate. However, the crust temperature would still remain high (§ 4.2) if the crust were 
very impure. Direct observational consequences of the core neutrino emissivity are unfortunately 
lacking. It is only in the cooling after accretion halts and the crust thermally relaxes (as in the 
transients; Brown, Bildsten, & Rutledge 1998) that the mode of core neutrino emissivity can be 
investigated. This is unlike the case of isolated, cooling neutron stars, for which the core neutrino 
cooling must be treated correctly. 
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The results of this investigation show that the most vexing impediment to further calculations 
of the thermal structure of an accreting neutron star, and hence to a better understanding of the 
issues raised in this section, is the need to calculate the composition throughout the crust for the 
trajectory in {n,T) space followed by an accreted fluid element. 

It is a pleasure to thank Lars Bildsten, Andrew Gumming, Andrew Melatos, and Greg Ushomirsky 
for many helpful discussions and for reading drafts of this work. I also thank Ghris Pethick for 
suggesting that the nuclei in the inner crust may remain spherical if the charge is low enough 
and the referee for helpful comments on the melting of a multi-species crystal. This research was 
supported by NASA grant NAG5-8658. EFB is supported by a NASA GSRP Graduate Fellowship 
under grant NGT5-50052. 
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